Fine Residual Carrier Frequency and Sampling 
Frequency Estimation in Wireless OFDM Systems 



Chen Chen, Yun Chen*, and Xiaoyang Zeng 
ASIC and System State Key Laboratory, Fudan University, Shanghai, China, 201203 

Email: chenyun@fudan.edu.cn 



(N 

o 

(N 

> 
O 

00 



<Z3 



> 
o 

CO 

oo 



CN 



X 



Abstract — This paper presents a novel algorithm for residual 
phase estimation in wireless OFDM systems, including the 
carrier frequency offset (CFO) and the sampling frequency 
offset (SFO). The subcarriers are partitioned into several regions 
which exhibit pairwise correlations. The phase increment between 
successive OFDM blocks is exploited which can be estimated 
by two estimators with different computational loads. Numerical 
results of estimation variance are presented. Simulations indicate 
performance improvement of the proposed technique over several 
conventional schemes in a multipath channel. 

I. Introduction 

Although the Orthogonal-Frequency-Division-Multiplexing 
(OFDM) technique significantly enhances the system perfor- 
mance under frequency-selective fading channels, it is vulner- 
able to synchronization non-idealities, including the symbol 
timing offset (STO), carrier frequency offset (CFO), and 
sampling frequency offset (SFO). 

The previous works including 0~|-||3] deal with the coarse 
STO and CFO estimation in time domain before Fast Fourier 
Transform (FFT). However, due to the imperfections of com- 
pensation, after FFT, the residual part of CFO remains to be 
corrected. Also, at this stage, SFO should be estimated and 
removed; otherwise, it would lead to a phase rotation not 
only proportional to the tone index within one OFDM block 
(inter-block increment), but also grows linearly for successive 
OFDM blocks (intra-block increment) J4] . 

In literature, several schemes are proposed to estimate or 
track the residual CFO and SFO in frequency domain with 
the assistance from pilot subcarriers |5)-0. In 0, Speth et 
al. utilize the symmetric locations of pilots to estimate CFO 
and SFO jointly. However, its performance degrades in the 
multipath channels. [6 | suggests three estimators with the help 
of the least square estimation (LSE). An improved weighted 
LSE variant is proposed by Tsai et al. in J7J which requires the 
second-order statistics of the channel state information (CSI). 
In general, these schemes mainly rely on the linearly growing 
inter-block increment. 

This paper proposes a novel technique to make use of the 
intra-block increment spanning a number of OFDM blocks. By 
dividing the subcarrier index into several regions, the method 
is capable of exploiting the pairwise correlation which leads to 
accurate results after applying least square fitting. Two variants 
differing in computation complexity are presented with their 
numerical variances derived. 

The rest of the paper is structured as follows. Section 
im presents the signal model in presence of CFO and SFO. 



Section |TTT] introduces the proposed technique and analytical 
results of variance. Simulation results are given in Section |IV] 
Finally, Section [V] concludes the paper. 

II. OFDM Signal Model with CFO and SFO 

We consider an OFDM system where the transmitted data 
is modulated by an iV-point Inverse FFT (IFFT). Assuming 
a total of M OFDM blocks to be delivered and each block 
consists of K data samples (K < N), the complex baseband 
signal is described by 
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where X are the locations of the K data subcarriers; for k ^ X, 
Xz fc is either pilot or null subcarrier. N g is the length of the 
guard interval, Nb the total length of an entire OFDM block 
given by Nb — N + N g , and T s the sampling interval, n(-) 
is the rectangular function defined as 



n(t) 



1 < t < N B T S 
otherwise 



The multipath channel is 

L-l 

h{t,r) = Y,ht(t)S(T-Te) 



(2) 



(3) 



where L is the total number of taps, {hi(t)}i=o,i,— the 
independent and Rayleigh distributed complex channel gains, 
{ri}e=o,i,--- the timing delay of each path, and 6(-) the 
delta function. Here, we assume that t# — IT S . 

Up-converting s(t) to a carrier frequency fx, the post- 
channel equivalent signal takes the form 



y{t) = [s{t)e^ fTt * h(t, r)] + w(t) 



(4) 



where the notation * stands for linear convolution, and w(t) 
the complex, identically independently distributed (i.i.d.), ad- 
ditive white Gaussian noise (AWGN) with zero mean and 
variance <r^; also, it is wide sense stationary (WSS), with 
independent real and imaginary part, and equal variance in 
both parts (er^/2). Now, assuming a CFO Af and a SFO r\ 
given as 



A/ ^ f T - f R 
rj 4 (T' s - T s )/T s 



(5) 
(6) 



where fn is the deviated carrier frequency and T' s the deviated 
sampling interval at the receiver. The received n-th sample in 
the Z-th OFDM block is 



n, 
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(7) 

After discarding the N g samples in the guard interval, the 
complex data for the Z-th block and on the fc-th subcarrier is 

m 
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where Hi , k is the channel transfer function (CTF) in frequency 
domain; 9^ ps e + rjk and e = AfNsT s the normalized CFO 
to the subcarrier spacing; a(Q k ) is the amplitude attenuation 
approaching unity and can be safely neglected; ICIj £ the 
inter-carrier interference (ICI) due to distorted orthogonality of 
subcarriers; Wi k is the WSS i.i.d. Gaussian noise in frequency 
domain with independent real and imaginary parts. Without 
loss of generality, e and r] can be regarded as the residual part 
of CFO and SFO after coarse synchronization or imperfect 
channel estimation and equalization. 

III. Proposed Technique 

Define the full set of subcarrier index asK = {fc|0<fc< 
N — 1}, which can be further divided into Q equally-spaced 
regions (assuming even N and Q, and N is divisible by Q), 
denoted as IK = K x U K 2 U • ■ ■ U K q U ■ • • U K Q where 
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(9) 



Here, Z denotes integers. Ignoring disturbances of ICI, using 
equation (O, for the q-th segment in the l-th OFDM block, 
the pairwise correlation is 

kl + k 2 =N q 
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and g — N g /N. Clearly, the extra phase rotation of the useful 
part in ( fT3l is irrelevant to subcarrier index k\ and ku", it is 
only pertinent to the OFDM block index I and segment index 



q. The cross terms are the main disturbance in estimation. In 
practice, the contribution of signal and channel (X l ki fc ) should 
be replaced by 



At 



where 
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The notation P is the full set of pilots and D the full set of 
null subcarriers. Xi k is the estimated X^ k , obtained by the 
decision feedback device. Hi k is the estimated CTF. V^ k fe2 
could be combined using a certain weight Ti :kl: k 2 given as 

J x [i^, fcl i 2 iH 1 ., 1 i 2 +ii. fe2 i 2 i« ! , fc2 i 2 J+^,,' Wei § hted 

ll, Simplified 

(18) 



where o\y = E{|W;,fc| 2 }. See Appendix [Al for the details 
of such selection for the weighted Ti k lt ka- Computation of 
weighted Ti^fa requires the second-order statistics of signal, 
channel and noise, avoided by the simplified scheme. For 
constant-modulus modulation, the weighted Ti t k lt k 2 reduces 
to 



?lMto = ([\ h im\ 2 + \Him\ 2 ] oWw+v 



(19) 



where er| = E{|X;^| 2 } = Const.. To obtain estimation of 



9 1 , we coherently stack V^ ki k2 by 
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where (■)* denotes the conjugation of its argument and 
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Kq.+ is the left half of K q while K 9 ._ the right one; U c is the 
absolute complement of U given by K\U = P U X; in general, 
the two-dimensional set C q ^ 0. B\ t can be estimated by 

arg{Zf}, q = l,2,--- ,Q 



The M x 1 vector 9 q = [6 q e v 9\ ^ 



linearized into 



■ • , Q (22) 
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where A is the M x 2 observation matrix expressed by 
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Fig. 1. The framework of the proposed technique. 8 denotes the extra phase 
rotation caused by CFO and SFO. Q = 2. The red dot and green dot form a 
pairwise correlation pair. 



The 2x1 vectors b q take the form 



fog — \* <i ' <: 



(27) 



where c q = 2e + r/N q ; x q lS the associated estimation error 
vector. By least square fitting, h q is given by 



h q = (A 1 A) 



A J 9 q (28) 



Note that both [b g ]i.i and [b 9 ]2,i give estimation of c q ; 

denotes the (i,j)-th entry of a vector/matrix. Here, we 
choose [b g ]i i and arrange all c q into the Q x 1 vector 
c = [ci C2 • • • c q ■ ■ ■ cq] t which leads to 

c = B f i + Xc (29) 
where fi — [rj e] T , x c * s me Q x 1 error vector, and 
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(30) 



is the Q x 2 observation matrix. Another least square fitting 
yields 

fi = (B T B) 1 B T c (31) 

The estimated rj and e" are and [/i]i,2 respectively. A 

simple sketch with Q = 2 is drawn in Fig. Q] Assuming 
correctness in tackling the phase ambiguity in the linearization 
process, derived in Appendix [A] estimation using either the 
weighted or simplified IY^ ^ is unbiased. On the other hand, 



the numerical variances of rj and e are 

Var {rj} 



81 Y^Q 77 (2Z-M+1) 2 

81 2^5=1^9 z^= — ft; — 



Var{e} 
where 

Fi, q = 

U a 



327V 4 (Q 2 - 1) 2 tt 2 (1 + g) 2 M 2 (M 2 
81 E? 



r^M-l (2i-Af + l) J 
2^i=0 F, „ 



327V 4 (Q 2 - 1) 2 tt 2 (1 + g) 2 M 2 (M 2 



1 



E(fcl,fc: 



(2Z(fc 1 ,k 2 )eC q r t'l,k 1 ,k2^ 



S(fc 1 ,i !2 ) e c, "^i.fo! ,fe 2 ^l,k 1 ,k 2 + 

*7fcl 



Weighted 
Simplified 



(32) 



(33) 



(34) 



(35) 



Lki 



IV 



l.k- 2 



1.M2 



w 



Y n = AN 



16V 2 [2g-Q-l] 2 
4Q 2 



2q - 1 



3Q 



Appendix |A] validates i*),^ Weighted > F^ q \ s 
Var {77} (weighted < 



imp. 



(36) 
(37) 
(38) 
and thus 



Var{r7}|si mp ., Var {e}| weighted 



< 



Var {e}|simp.- The equality establishes if and only if (iff) 
Al: the channel experiences flat fading {\Hi t k\ 2 = 1/iV) 
A2: constant modulus modulation ( | , | 2 = f|) 
Otherwise, weighted estimation always outperforms simplified 
estimation. To achieve the best performance, further assuming 
A3: equal and maximal cardinality of each set C q , de- 
noted by N Cq = N/2Q 
The variances under conditions A1^A3 are 

18Q 2 



Var {rj} 



Var {7} = 



tt 2 (1 + g) 2 M(M + 1)(M - l)iV 3 (Q 2 - 1)SNR 

(39) 

6(4Q 2 - 1) 

4tt 2 (1 + g) 2 M(M + 1)(M - 1)N(Q 2 - 1)SNR 

(40) 



where SNR 4 a 2 a 2 H /a 2 w ; a 2 H = E{\H^ k \ 2 }. 
Remarks: 

i) Apparently, an immediate way to enhance the perfor- 
mance is to raise M, which leads to asymptotically decreasing 
variances in cubic scale. Nevertheless, if the application is 
real-time oriented rather than quality preferred, where r\ and 
e should be tracked in the fastest manner, M and Q should be 
replaced by their minimums as min{Af} = 2, min{Q} = 2. 

ii) If [b g ]2.i is used for estimation, (|39l and (l40b are 
rewritten into 

Var {rj}' = W ~ *)<M ~W + 9?^ { ~ } (41) 



(1 + 2g) 2 

Var { er=^ M -y M -y i + ^Var { ^} 



(42) 



(1 + 2 ff ) 2 

which are significantly higher than d39b and d40b increasing 
squarely with M. Therefore, it is reasonable to use [b g ]i i. 



iii) According to ( fTOb , statistically, the proposed technique 
only relies on the independency and equal variance of the real 
and imaginary parts, and WSS assumptions of noise; it does 
not require the power spectrum density (PSD) of noise to be 
strictly fiat (white), since the expectation of the cross terms is 
zero. 

IV. Simulation 

In this section, we consider a wireless OFDM system with 
FFT size N = 512. The guard interval is N g = 64. Thus, 
the length of an entire OFDM block is N B = 576. The total 
number of OFDM blocks is M = 10, and the total number 
of segments is Q — 4. The carrier frequency is set at 5 
GHz. The sampling period is T s = 100ns. For brevity and to 
exploit the best performance of the proposed estimators as well 
as other conventional pilot-assisted schemes, all subcarriers 
are regarded as pilots; otherwise, notations in (fTTI i must be 
used which varies with the accuracy of decision feedback 
device. The signal is modulated from 16-PSK constellation. 
The channel consists of L = 32 Rayleigh taps, which are 
statistically independent distributed with a power delay profile 
decaying exponentially: 

exp(-£/L) 



Et"o |ex P H/i)| s 



£ = 0,1,2,- ■■ ,L- 1 



(43) 



For the proposed estimators and the scheme in J7J, CTF 
is assumed to be known perfectly as well as cr^/LJ. unless 
otherwise mentioned. Mean squared error (MSE) results are 
used to benchmark the performance, defined as MSEj?}} = 
E {\V-V\ 2 } and MSE{?} = E{|?-e| 2 } where E{-} denotes 
the expectation of its argument. 

A. Comparison o/MSE{r/} 

Fig. |2] highlights the comparison of MSE{f/} among the 
proposed estimators with other schemes in f5]-|l7). Numerical 
result of Varjry} in (l39l is drawn by assuming A1^A3. In 
the multipath channel, both of the weighted and simplified 
estimator achieve the best performances. In the flat fading 
scenario, the weighted estimator reduces to the simplified one. 
Var {rj} provides a tight bound in moderate SNR. 

B. Comparison o/MSE{e} 

Fig. [3] shows the performance comparison of MSE{e} in 
multipath channel. [5| could achieve the best performance 
when SNR < 7 dB, which cannot be sustained into higher 
SNR. Again, Var {e} provides a tight bound in moderate SNR. 

C. MSE{ry} with a Varying rj 

Fig. |4] shows the deviation of MSEjry} when rj changes 
under SNR = 20 dB in multipath channel. For the proposed 
estimators, MSE{r/} is asymmetric for negative and positive 
rj due to the presence of a positive e; for rj > 0, the 
performance degrades gradually with a higher rj since the ICI 
is increasing simultaneously. For a major part, |7| and the 
simplified estimator entangle with each other. 

'For OFDM systems containing null subcarriers, <xjL could be estimated, 
which is omitted in this paper. 



D. MSE{e} with a Varying e 

Fig. [5] displays the deviation of MSE{e} with a varying 
e under SNR = 20 dB in multipath channel. Different 
from Fig. |4j shape of MSE{e} is akin to symmetric, since 
comparing with e, the contribution of rj on ICI is minor. In the 
full range, both of the proposed schemes outperform others. 

E. MSE{^} with Different Estimation Accuracy 

Fig. [6] shows the MSE{?7} performance with different 
channel estimation accuracy in multipath channel. Employing 
the same philosophy of (9 J, the inaccurate channel estimation 
takes the form 



Vl — K 2 Him + KJlk 



(44) 



where k represents the estimation accuracy, and Ji k the 
additional complex noise with zero mean and variance 0.5 
in its real and imaginary part respectively, independent from 
H^k- The weighted estimator degrades significantly when 
severe inaccuracy occurs, which hinders the performance 
improvements especially in moderate to high SNR region. For 
moderate to high SNR, the weighted estimator outperforms 
the simplified one under different k. Similar conclusion can 
be drawn for MSE{?}. 

F. MSE{?7} and MSE{e} under Mobility 

Fig. exhibits the MSE{?7} and MSE{?} in presence 
of terminal mobility with SNR = 20 dB under multipath 
channel. Merely the CSI pertinent to the first moment of the 
Rayleigh fading channel in simulation (t=Q) is assumed to 
be known a priori; for the ensuing frames, the same CSI 
is used which entails a loss in channel estimation accuracy. 
The Doppler bandwidths with respect to terminal speed of 
[50, 100, 150, 200]km/h are 232, 463, 695, 927 Hz. In general, 
the performance deviation is insignificant if not imperceptible 
even the terminal velocity reaches 200km/ h, since the maxi- 
mal value of the product between the Doppler bandwidth and 
the duration of an OFDM block is 5.3 x 10~ 2 , a relatively small 
value. Thus, the CSI is sound enough to secure an excellent 
estimation. 

V. Concluding Remarks 

In this paper, we propose a joint estimation technique to 
deal with residual CFO and SFO estimation. By dividing the 
subcarrier index into a number of regions and exploiting the 
pairwise correlation, we estimate the phase increment between 
adjacent OFDM blocks, which yields accurate estimation after 
Q + 1 times of least square fitting. Extensive simulations 
indicate better performance over several conventional pilot- 
assisted schemes. 

Appendix A 
Bias and Variance of Estimation 

First of all, consider the case of the simplified Tj^,^ in 
([Tol l and the statistical information of 6 q in ( f22b . For the Z-th 
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Fig. 4. MSE{f;} under different r/; M = 10, Q = 4, r) G [-10 : 10] X Fig. 5. MSE{?} under different e; M = 10, Q = 4, T? = 5 X 10 _5 ,e S 
10 _5 ,e = 0.02, and SNR = 20 dB. [-0.1 : 0.1], and SNR = 20 dB. 



element of q {Of e ), subtracting the left hand side of ( f22b 
yields 



A# 9 - tf 9 



^ = arg[Zf/exp(j^^)] (A.45) 



If 0f e _ is in vicinity of # 9 e , using the approximation 
tan(:r) ps a; for a; small enough, we may write 



c itA8 q \ 
tan(A^ e v ) « ^ , ^ ~ At? 9 



l,e,r/ 



(A.46) 



where 3t{-} and 5{ } represent the real and imaginary part of 
the arguments. Expectation of Af? 9 e in (IA.46I ) is 



E 



/ S{Ac^} | E{S{A0f [6[ „}} 



which holds if E{3?[A6» 9 et; ]} > [Var^Afl?^]}] 1 / 2 . 
In fact, each component in Sy{A# g £ } contains either 
Wi M , W iM , or W lM W lM (see and therefore, 

E{$j{A0? }} = 0, which finally leads to the unbiasedness 



of fi in ( Bil l since only linear intermediate operations are 
involved. 



For the numerical variance of A# 9 e , we could use 



Var [A6f 



V a r[Z{A0l n }} 



(A.48) 



if E{3ft[Aef e ^]} > [Vax^Aflf^]}] 1 / 2 . Standard calcula- 
tions yield 



Var[3{A0 9 e „}] = £ a 2 w \H lM \ 2 \H l<k2 f 

ki ,k 2 

(A.47) [\X lM \ 2 \H lM \ 2 + \X lM \ 2 \H lM \ 2 + a 2 v ] 



(A.49) 
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Fig. 6. MSE{»7} under different accuracy of channel estimation; k 6 [0.1 : Fig. 7. MSEjr;} and MSE{e} under different terminal velocity; M = 
0.5], M = 10, Q = 4,r; = 5 X 10" 5 ,e = 0.02, and SNR 6 [0,30] dB. 10, Q = 4,?? = 5 X 10" 5 ,e = 0.02, and SNR = 20 dB. 



(E[5ft{A^ er; }]) 2 = Yl \X lM \ 2 \X lM \ 2 \H lM \*\H lM 

^ ki,k 2 

(A.50) 

Note that, for visual clearance, we abbreviate the notation 
(ki, k 2 ) £ C, with ki,k 2 . Thus, 



Var[Aflf ( J«fi,,| slmp . 



(A.51) 



where Fi q is defined in ( [34) . Using the linear intermediate 
manipulations, we derive (l32b and Q31 >. Substituting i^^lsimp. 
with i 7 /^ | weighted, backward produces the weighted version of 
^i,k ± ,k 2 m < fl~8T > after lengthy calculations. 

To prove F /:g | Woightc d > F; i9 | Simp ., we invoke the Cauchy- 
Schwarz-Inequality [ 10 1 . The essential steps are listed below. 
Proof: To prove 



E 



/,fel,fe 2 ^-i—ik\,k 2 ^l,k\,k 2 > (A 52) 



is equivalent to prove 



(A.53) 



Now, using i4(fci,*a) = and B(fci,*a) = 

E fcl , fc2 ^IMM^MM + X )' the Cauchy-Schwarz-Inequality 
gives 



and we have the conditions Al, A2 0. Therefore, we verify 

that F^g (weighted > -F/,qr|simp.- ■ 
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Mh,k 2 ) Y B(h,k 2 ) > ! Y y/A(k 1 ,k 2 )B(k 1 ,k 2 ) 
fci,fc 2 fci,fc 2 I fci,fc 2 

(A.54) 

where the right hand side of (IA.54) is exactly the right hand 
side of (TA331 . The '=' holds iff 

tlMte = Const '' V ( fc i> fc 2) € C q (A.55) 



2 Actually, it is merely a sufficient condition of jA.55l . Specifically, consider 
that there is a sage at the transmitter side who could render \Xi^\ 2 = 
G\Hi fe| — 2 , Vfc where C is a constant, then jA.55l also establishes. However, 
it is not insightful to be pursued. 



